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Abstract 

This is a pedagogical introduction to the basics of the cosmic string 
theory and also a review of the recent progress made with respect 
to the macrophysical predictions of the theory. Topics covered in- 
clude, string formation and evolution, large scale structure formation, 
generation of peculiar velocity flows, cosmic microwave background 
(CMB) fluctuations, lensing, gravitational waves and constraints from 
the msec pulsar. Particular emphasis is placed on the signatures pre- 
dicted on the CMB and the corresponding non-gaussian features. 
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1 Introduction-Cosmology Basics 



The Big Bang theory|jT[| is the best theory we have for describing the universe. 
It is a particularly simple and highly successful theory. It makes three basic 
observationally consistent assumptions and derives from them some highly 
non-trivial predictions which have been verified by observations. The main 
assumptions are the following: 

• The universe is a homogeneous and isotropic thermal bath of the known 
particles (cosmological principle 0). 

• General Relativity is the correct theory to describe physics on cosmo- 
logical scales. 

• The energy momentum tensor of the universe is well approximated by 
that of a perfect fluid i.e. 

Tij = diag[p, -p, -p, -p] (1) 



The first two assumptions imply that the universe can be described by the 
Robert son- Walker metric of the form 

ds^ = de - a{tf[ ^ + r\d9^ - sin^ edcf)^)] (2) 

where a{t) is known as the scale factor of the universe and k can be normal- 
ized to the values —1, 0, 1 for an open (negative curvature), spatially flat and 
closed (positive curvature) universe respectively. 

The third assumption can be used along with the Einstein equations to 
derive the Friedman equations that determine the dynamics of the scale factor 
a{t). These may be written as 

(2 2 k 8nG 

y - ^=^^ 

a AnG , , , . 

- = —iP + 3p) (4) 

These are two equations with three unknown functions (a{t) , p(t) , p{t)) 
and therefore the equation of state p = p{p) is also needed in order to obtain 
a solution. For example, for a flat [k = 0) radiation dominated universe 
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(p = P/3), it is easy to show that a(t) ~ t^^"^ while for a matter dominated 
universe {p = 0) we have a{t) ~ t^/^. 

Observations indicate that at the present time the matter density domi- 
nates over radiation density in the universe {pmatito) ^ Pradito)) and there- 
fore we hve in a matter dominated universe. However, this was not the case 
at early times. The Friedman equations may be used to show that 



i(pa') = -3aV (5) 



which implies 



P.^atit) ~ a(t)-3 (6) 
Prad{t) ~ a{t)-^ (7) 

and therefore there is a time t^q < to of equal matter and radiation such that 

t <teq^ Pmat{t) < Prad{t) (8) 

Using Eq. (7) and the well known result from statistical mechanics that 
Prad ~ T'^ we obtain that the temperature of radiation scales as 

T(t)~a(t)-i (9) 

in a matter or radiation dominated universe. Cosmic Microwave Background 
(CMB) measurements^ have shown that T(to) = 2.736±0.017 and therefore 
the temperature at earlier times is predicted to be 

Tit) = 2.7'>K^ (10) 

The three main predictions of the Big Bang theory are the following: 

• The expansion of the universe (Bubble's law). 

• The existence of the Cosmic Microwave Background (CMB). 

• The relative abundances of the light elements (nucleosynthesis). 
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A simple way to derive Hubble's law in the case of a flat (k=0) universe is 
to consider two points in space separated by coordinate (comoving) distance 
dr (with d6 = d(j) = 0). The proper (physical) distance between the two 
points is given by 

ds = a{t)dr (11) 
Therefore, the physical relative velocity between the two points is 

s = dr + af (12) 

where dr is known as the Hubble flow and ar is the peculiar velocity. For 
large comoving distances r the second term can be ignored and we have 

f = s ~ -s = H{t)s ^ ~ w = H{t)s (13) 
ct 

where z is the Doppler redshift and H{t) = ^ is the Hubble constant at 
cosmic time t. Eq. (13) is known as the Hubble's law^ and has been verified 
by observations. The value of -ff(to) may be written as 

H{to) = lOO/i km/{sec ■ Mpc) (14) 

with he [1/2,1]. 

A particularly important prediction of the Big Bang theory is the ex- 
istence of the CMB|^. Penzias and Wilson observed for the first time 
in 1965 a highly isotropic backgroung of microwave photons with spectrum 
that was thermal to a high accuracy. Several experiments since then have 
verified this observation and today it is known that the temperature of this 
background is given by Eq. (10) {t = to) and temperature anisotropics are of 
O(10~^). An additional term of dipole anisotropy ||^ is known to be present 
due to our motion with respect to the CMB. The magnitude of this term is 
^ ~ 10~^ and corresponds to a velocity of the earth of about 600fcm/sec 
with respect to the CMB frame. Even though the dipole term dominates 
over the primordial fluctuations it can be easily subtracted due to its dipole 
nature. 

CMB photons are free at the present time to (they have mean free path 
of cosmological scale) and there is no apparent mechanism that could have 
caused their thermalization. However, such mechanism is naturally pro- 
vided by the Big Bang theory. According to this theory, the universe is 
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hotter at early times (cf Eq. (10)) and therefore there is a time tree < to 
{time of recombination) when photons are energetic enough to ionize matter 
{T{tree) ~ ^ion — O(IO^) ^ K) . In such a primordial plasma, thermalization 
of photons can occur effect ively[|], ||. 

The redshift z{tree) at the time of recombination defined as 

1 , N _ _a{to) Tjtree) (^^s 

' + ^<'"'> = ^ = W 

is therefore z{tree) — 1500. Thus, according to the Big Bang theory the CMB 
photons were scattered for last time at t = tree and carry a 'photograph' of 
the universe taken when it was 1500 times smaller and hotter. There is 
exciting information hidden in this map. An as yet unknown source has 
created primordial fluctuations that evolved gravitationally into what we see 
as galaxies, clusters and large scale structure today. It is those fluctuations 
in their primordial gravitationally unaffected form that have been recorded 
by the CMB photons. Clearly, information on such a primordial pattern 
can impose severe constraints on theories attempting to explain the origin of 
primordial fluctuations. 

The third prediction of the Big Bang theory is made on the relative 
cosmological abundances of the light elements {H, He, H^, Li"^) which span 
10 orders of magnitude and are predicted correctly by the Big Bang theory 
provided that the baryon to photon ratio is ?7 = ^ ~ 10~^°|T0|. 

The assumptions of homogeneity and isotropy made by the Big Bang the- 
ory, even though consistent with very large scale (~ 100h~^Mpc) observations |TT 
are only first approximations to the realistic system. It is enough to look at 
the night sky to realize that the universe is not homogeneous and isotropic on 
relatively small scales. Instead, matter tends to cluster gravitationally and 
form well defined structures. In fact, detailed large scale structure observa- 
tions have indicated the existence of sheets and filaments of galaxies sepa- 
rated by large voids. The typical scale of these structures is about 50h~^Mpc 



i2,ra,ra, 15 



Fig. 1 shows a slice of the Cf A survey which is a map of galaxies in redshift 
space (depth in the sky). The presence of the above described structures is 
evident in this map (the axis of redshifts can easily be converted to distance 
(depth in the sky) by using Bubble's law|^). Can these large structures be 

■^For example for h ^ 1 dividing z by 100 gives the corresponding distance in Mpc. 
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explained by assuming that matter moved due to gravity from an initially 
homogeneous state? 

Figure 1: The CfA Survey shows a map of galaxies in redshift space 
(depth in the sky). This 'slice of the universe' includes a range of about 6° 
in declination space. 



The maximum distance that matter can have travelled since the Big Bang 

is 

Srmax ~ ^i; to - 6h~^Mpc (16) 
where Sv ^ 600km / sec is the typical peculiar velocity of galaxies (induced by 
gravity) at the present time to- This dr,nax is much less than the typical scale 
of observed structures. Thus, these structures can not be the result of gravity 
alone. They must also reflect the presence of primordial perturbations. The 
question that we want to address is What caused these fluctuations? 

Until about 15 years ago there was no physically motivated mechanism to 
cause these perturbations. However, during the past decade, two classes of 
theories motivated from quantum field theory have emerged and attempt to 
give physically motivated answers to the question of the origin of structure 
in the universe. 

A model for large scale structure formation is characterized by two basic 
features. The first is the kind of the assumed dark matter and the second is 
the type of primordial fluctuations. 
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Dynamical measurements based mainly on peculiar velocities have shown 



that at least 90% of the matter in the universe is not luminus|16, 17, 18|. This 
non-luminus matter which has not been directly detected yet, is called dark 
matter and its properties determine critically the rate of gravitational growth 
of primordial perturbations at different cosmic epochs. In particular, dark 
matter with low non-relativistic velocities at the time teq (when fluctuations 



start to grow) is called cold dark matter (CDM)[0. A CDM candidate 



motivated from particle physics is the hypothetical particle axzon||20[ which 
appears to be a necessary consequence of a successful resolution of the strong 
CP problem in QCD. Dark matter which has relativistic velocities at tgq is 
called hot dark matter (HDM) and a primary candidate for it is a massive 
{25h~'^eV) neutrino. Due to their relativistic velocities HDM particles can 
not cluster gravitationally at early times and small scales. This effect is called 
free streaming^^. The critical scale for HDM perturbation growth is the free 



streaming scale lfs{t) ^ v^{t)t i.e. the distance travelled by a HDM particle 
in a Hubble time t. Adiabatic perturbations (those produced by inflation) 
on scales I <lfs are erased on small (galactic) scales due to the effects of free 
streaming. 

The second characteristic of large scale structure formation models is the 
type of primordial fluctuations. There are two broad classes of primordial 
perturbations, both produced by physically motivated mechanisms. The first 
includes gaussian adiabatic perturbations produced ||22|, ^ |25[ during a 



period of exponential growth of the universe known as inflation^^. These 
perturbations may be represented as a superposition of plane waves with 
random phases i. e. 

5 = ^(x) ~^|4|e*''e*'^" (17) 
P k 

where ^(x) represents the primordial density fluctuation pattern and 9k 
are random uncorrelated phases. By the central limit theorem, since 5 is 
a superposition of an infinite number of uncorrelated random variables, its 
probability distribution is gaussian i.e. 

P{6) ~ e-'^'/'^' (18) 

Fig. 2a shows schematically the features of such a gaussian perturbation 
pattern. 
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Figure 2: Gaussian fluctuations may be represented as a superposition 
of plane waves with random phases (a) while topological defect fluctuations 
are represented by a superposition of seed functions (b). 



The second class includes primordial perturbations produced by a su- 
perposition of localized seeds created during a phase transition in the early 
universe. Such seeds are known as topological defects^T^ ^ ^ ^ ^ and 



their formation is predicted by many (but not all) Grand Unified Theories 
(GUTs)[0. The pattern of topological defect perturbations may be repre- 
sented as a superposition of localized functions (Fig. 2b) i.e. 

6 = J:^Ax-x.) (19) 

i P 

and therefore, in general the probability distribution of S is not gaussian^) 

Both inflationary and topological defect perturbations can be combined 
with either HDM or CDM for the construction of structure formation theo- 
ries. In models based on gaussian adiabatic fluctuations with HDM, galactic 
scale fluctuations are erased by free streaming and therefore these structures 
can only form at later times by fragmentation of larger objects. This is in 
contrast with observations indicating that galaxies formed earlier than larger 



*For a not very large number of superposed seeds. 
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structures. These models face also other conflicts with observations {eg they 
violate CMB fluctuation constraints) and have therefore been placed into 
disfavor. On the other hand, the model based on inflationary gaussian fluc- 
tuations with CDM has been well studied and makes concrete predictions 
that are in reasonable agreement with most types of observations especially 
on small and intermediate scales. This model is currently the 'standard 



model' for structure formation 33 



Models based on topological defect perturbations combined with either 
hot or cold dark matter have not been as well studied but they have been 
shown to have several interesting features that make them worth of further 
investigation. Most of the remaining of these talks will focus on one of the 
most interesting types of topological defects cosmic strings^^, Three 
basic aspects of cosmic string physics will be discussed: their formation, 
evolution and gravitational effects. 



2 Formation of Topological Defects 
2.1 Kibble Mechanism 

Even though no particle corresponding to a scalar field has been discovered 
so far, the role of scalar fields in particle physics models is central. Indeed, 
scalar fields provide a natural and simple mechanism to induce spontaneous 
symmetry breaking in gauge theories thus achieving two important goals: 
First maintain gauge invariance and renormalizability of these theories and 
second give mass to gauge bosons thus making them consistent with inducing 
short range interactions like the strong and electroweak forces. Symmetry 
breaking can be achieved by the use of scalar field potentials of the form 

= A($t$ _^2)2 (20) 

where $ = [$1, $Ar] is a multiplet of scalar fields. 
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Figure 3: The Kibble Mechanism: The evolution of a scalar field in a 
potential with discrete minima ((a) left) leads to the formation of a domain 
wall in physical space after the relaxation of the field in its minima ((a) 
right). When the potential minimum (vacuum manifold) has the topology 
of a circle (S^) a cosmic string forms in physical space (b) while a monopole 
forms when the vacuum manifold is a sphere (5*^) (c). 



Consider for example a one component scalar field whose dynamics is 
determined by the double well potential of Fig. 3a. At early times and high 
temperature T the field $ will have enough energy to span the whole range 
of the potential^ and go over the barrier between the two minima. As the 
universe cools and expands however, the energy of the scalar field will drop 

^The temperature dependence of the effective potential is ignored here since its calcu- 
lation is highly non-trivial in out of equilibrium systems like the one discussed here. 
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and it will eventually be confined in one of the two minima [+1], —rj). The 
choice of the minimum is arbitrary. In fact, there will be neighbouring un- 
correlated or causally disconnected regions where the choice of the minimum 
(vacuum) will be opposite (Fig. 3a). By continuity of the field $ there will 
be a surface separating the two domains where $ = and therefore $ will be 
associated with the high energy of the local potential maximum. This surface 
of trapped energy density is known as a domain wall and is the simplest type 
of topological defect. Domain walls form in theories where the vacuum mani- 
fold M (minimum of potential) has more than one disconnected components. 
This fact is expressed in homotopy theory [p6| as tiq^M) ^ 1. 

The mechanism described above for the formation of domain walls is 
known as the Kibble mechanismETi 123] and applies in a similar way to the 



case of other topological defects. Consider for example the case where the 
dynamics of a two component scalar is determined by the potential of Fig. 
3b. After the field relaxes to its vacuum, there will be (by causality) regions 
of the universe where the field $ will span the whole vacuum manifold as we 
travel around a circle in physical space. Topologically stable vortices can also 
form in systems with multiple scalars[^ while metastable vortices can even 
form in systems where the vacuum manifold is not S^. In the later cases the 
stability is not topological by dynamical ||38||. By continuity of $ there will be 
a point inside this circle where $ = 0. Such a point (and its neighbourhood) 
will be associated with high energy density due to the local maximum of the 
potential at $ = 0. By extending this argument to three dimensions, the 
point becomes a line of trapped energy density known as the cosmic string. 
In general cosmic strings form in field theories where there are closed loops 
in the vacuum manifold M that can not be shrank to a point without leaving 
M. In homotopy theory terms, we require that the first homotopy group of 
the vacuum should be non-trivial i.e. 7ri(M) 7^ 1. It may be shown that this 
is equivalent to ttqIH) 7^ 1 where H is the unbroken group in the symmetry 
breaking G —>■ H (M=G/II). As will be discussed below, the number of 
times the field $ winds around the vacuum manifold as we span a circle in 
physical space around the string is a topologically conserved quantity called 
the topological charge (or winding number) and its conservation guarantees 
the stability of the string. 

In theories where the vacuum manifold is a sphere, the defect that forms 
by the Kibble mechanism is pointlike and is known as the monopole (Fig. 
3c). In general monopoles form in theories where the vacuum manifold has 
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unshrinkable spheres i.e. 7r2(M) 7^ (equivalent to vri(if) 7^ 1). Monopoles, 
like strings have associated topological charge whose conservation guarantees 
their stability. As will be discussed later, most types of domain walls and 
monopoles are inconsistent with the standard Big Bang theory as they lead 
to quick domination of the universe by defect energy and subsequent early 
collapse. On the other hand, in the case of strings there is a natural mech- 
anism to effectively convert string energy to radiation and thus prevent the 
string network from dominating the universe. 

Topological defects similar to the ones described above can form not only 
in a cosmological setup but also in several condensed matter systems like 
liquid crystals ICT, He^ |^ and superconductors!^ . 



2.2 Non-local defects: Textures 

All the defects discussed so far have two important properties: 

• The field configuration can not be continously deformed to the trivial 
zero energy vacuum where the field would point to the same direction 
everywhere in space. 

• There is a well defined region in space where most of the defect energy 
is localized. It will be shown that the size of this domain is proportional 
to the inverse of the symmetry breaking scale 77. 

From the second property it follows that the above described defects are 
localized defects. There is a second class of defects that have the first but not 



the second of the above properties. Textures belong to this later class 1 43, 44 
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Figure 4: The field configuration of a non-local texture in one (a) two 
(b) and three (c) space dimensions. 



Consider a global (no gauge fields) theory where the vacuum manifold is 
(a circle). Consider also the one dimensional (fixed boundary conditions) 
scalar field configuration shown in Fig. 4a. This configuration is topological 
(can not be deformed to the trivial vacuum) but is non-local (there is no core 
of energy density). It is known as the one dimensional texture and in general 
it is non-static. It is straightforward to generalize the one dimensional texture 
to two dimensions and three dimensions by considering theories with M = S"^ 



and M = S (Fig. 4b and 4c). In fact it may be shown (Derrick's theorem 



that a rescaling of the radial coordinate r by a scale factor a > 1 
(r — > ar) leads to a rescaling of the energy of the three dimensional texture by 
the inverse factor a {E E/a) which implies that collapse and subsequent 
unwinding is an energetically favored process provided that the topological 
charge is larger than a critical valuefl^, ^ Even though the physics 



and observational effects of textures are particularly interesting|51, 52, 53 



subjects they are outside of the scope of this review and will not be further 
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discussed here (for a good review of textures see Ref. p5|). 



2.3 Classical Field Theory 

In what follows I will focus on the cosmological effects of cosmic strings. The 
simplest model in which strings involving gauge fileds can form||5^ is the 
Abelian Higgs model involving the breaking of a U{1) gauge symmetry. The 
Lagrangian of the Abelian Higgs model is of the form 

C = ^D^$D^<I> - V{^) + \f^uF^'' (21) 

where F^, = d^A, - d,A^, D^ = d^- leA^ and V{^) = f (|$|2 - ri^)\ 
($ = $1 + i$2)- The field configuration corresponding to a stringjQ ^ 
in this theory is shown in Fig. 3b. The topological stability of the string 
(vortex in two space dimensions) is due to the conservation of the topological 
invariant (topological charge or winding number) 

da , . 

•" = /„ Tw"' (22' 

where a is an angular variable determining the orientation of $ in the vac- 
uum manifold which is (circle) in the Abelian Higgs model and 6 is the 
azimouthal angle around the vortex in physical space. For the string shown 
in Fig. 3b we have m = 1 and the field winds once in the positive direction 
as we span a circle around the string. 

The field configuration of a vortex may be described by the following 
ansatz: 

$ = /(r)e^'^^ (23) 

A = ^e. (24) 
r 

with /(O) = f(0) = 0, /(oo) = 7] and f(oo) = m/e. The forms of /(r) and 
v{r) can be obtained numerically ||56|| from the field equations with the above 



ansatz. The energy of the vortex configuration (energy of string per unit 
length) is 

, ^ f ^fj'.ir- + + i^^f' + _ ,ff) (26) 
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Clearly, t> 7^ is required for finite energy i.e. the gauge field is needed in 
order to screen the logarithmically divergent energy coming from the angular 
gradient of the scalar field. A vortex with no gauge fields is known as a global 
vortex and its energy per unit length diverges logarithmically with distance 
from the string core. This divergence however is not necessarily a problem in 
systems where there is a built-in scale cutoff like systems involving interacting 
vortices. In such systems the intervortex distance provides a natural cutoff 
scale for the energy integral. Similar considerations apply to monopoles with 
no gauge fields even though the energy divergence in that case is linear rather 
than logarithmic p7|, p^ , |59[] 

Even though no analytic solution has been found for the functions /(r) 
and v{r) it is straightforward to obtain the asymptotic form of these functions 
from the field equations (Nielsen-Olesen equations). The obtained asymp- 
totic form for r ^ 00 is^, ^ (but see Ref. for a correction to the 
standard result) 

/(r) ^ r^-^e-^'"^ (26) 

171 

v{r) — -c^v^e""'"' (27) 

where Cy and c/ are constants. Therefore the width of the vortex is 

w ~ T]-^ (28) 

where rj is the parameter of the symmetry breaking potential also known as 
the scale of symmetry breaking. The energy per unit length /i may also be 
approximated in terms of rj as 

/i ~ /" d'^xV{0) ~ 7]^ (29) 



The typical symmetry breaking scale for GUTs is 77 ~ lO^^GeV which leads 
to extremely thin and massive strings 

w ~ 10"^°cm (30) 
fi ~ lO^Hons/mm (31) 

Therefore, the Nielsen-Olesen vortex may be viewed as a linear, topo- 
logically stable field configuration involving both magnetic field energy and 
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Higgs potential energy that decay exponentially at large distances. No fi- 
nite gauge transformation can gauge away a vortex gauge field everywhere 
in space. The several interesting quantum field theoretical interactions of 
vortices with particles (fermions) are out of the scope of this review [|61| . 



3 Cosmic String Evolution 
3.1 Nambu Action 

The exact analytic treatment of the evolution of a network of strings would 
involve the analytic solution of the time dependent nonlinear field equations 
derived from the Lagrangian of Eq. (21) with arbitrary initial conditions. 
Unfortunately no analytic solution is known for these equations even for the 
simplest nontrivial ansatz of the static Nielsen-Olesen vortex. Even the nu- 
merical solution of these equations is practically impossible for systems of 
cosmological scales and with complicated initial conditions. The obvious al- 
ternative is to resort to realistic approximations in the numerical solution of 
the field equations in cosmological systems. First, the correct initial condi- 
tions for a numerical simulation must be obtained by simulating the above 
described Kibble mechanism on a lattice. This may be achieved by im- 
plementing a Monte-Carlo simulation implemented first by Vachaspati and 
Vilenkin||62||. They considered a discretization of the vacuum manifold (S^) 
into three points (Fig. 5) and then assigned randomly these three discrete 
phases to points on a square lattice in physical space in three dimensions. 
For those square plaquettes for which a complete winding of the phase in the 
vacuum manifold occured, they assigned a string segment passing through. 
It may be shown that this algorithm leads to strings with no ends within 
the lattice volume i.e. strings either form loops or go through the entire lat- 
tice volume (infinite strings). This simulation showed that the initial string 
network consists of 80% long strings and 20% loops. 
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Figure 5: The Vachaspati-Vilenkin algorithm used to study the forma- 
tion of a string network by Monte-Carlo simulations. 



In order to find the cosmological effects of this initial string network, it 
must be evolved in time. Since it is impractical to evolve the full field equa- 
tions on the vast range of scales 10~^° cm (string width) to lO^Mpc (largest 
cosmological scales) we must resort to some approximation scheme. The ratio 
w/R of the width of the string over the string coherence scale is an intrinsi- 
cally small parameter for cosmological strings {w/R << 10~^^ for t > teq)- 



Thus, this parameter can be used to develop a perturbation expansion [pSj |64| 
for the action that describes the dynamics of a string segment. The zeroth 
order approximation is obtained for w/R = 0. 

A heuristic derivation of the action of a zero width string may be obtained 
by using the analogy with a relativistic free point particle with mass m and 
velocity v. The relativistic action describing the motion of such a particle is 
of the form 

pB j-B I 1 

S = -m dT = -m dt{l - f 2)^/2 ~ / dt{-mv^) + const (32) 

J A J A J A 2 

where A, B are the endpoints of the trajectory and r is the particle's proper 
time. As shown in Eq. (32), the non- relativistic limit is obtained for low 
velocities. 
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The trajectory of the string is not parametrized by only one variable (the 
proper time r) but by two variables (a timelike variable r and a spacelike 
variable / which may be viewed as parametrizing the length of the string). 
Thus, the generalization of the point particle relativistic action to the string 
with mass per unit length /i is the Namhu action defined as: 

S = -fxJ^ dldr^f^ + 0{w/R) (33) 

where g'^'^'> = det{X^^X\g^y) (a,b = 1,2 —» (r, /), /i, = 1,...4) is the metric 
of the world sheet spanned by the string. The Nambu action is an excellent 
approximation to the dynamics of non-intersecting cosmic string segments 
and is much simpler to handle numerically than the full field theoretic action. 



Figure 6: The intercommutation (exchange of string segments) occurs 
when string segments intersect and favors formation of string loops (a). The 
effect is closely related to the right angle scattering of vortices in head on 
collisions (b,c) (from Ref. p7|). 

A crucial assumption in the derivation of the Nambu action is that string 
segments do not interact with each other. This assumption breaks down when 
two string segments intersect. The outcome of such an event can only be 
found by evolving the full field equations. Such numerical experiments have 
shown p5|, |66| , |67| that at intersections string segments exchange parteners 
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and a process known as intercommuting occurs (Fig. 6a). Intercommuting is 
closely related to the right angle scattering of vortices in head on collisions 
(Fig. 6b, c) which has been observed in several numerical experiments and can 
be analytically understood using dynamics in moduli spaces|6^. As shown 



in Fig. 6a intercommuting tends to favor the formation of string loops which 
oscillate and decay by emitting gravitational radiation with characteristic 
frequency u ~ {R is the loop radius) and rate E = •jGfi'^ (7 = const ~ 
50)11. 



Therefore, this mechanism for loop formation, provides also an efficient 
way for converting long string energy that redshifts with the universe expan- 
sion as Estr ~ to radiation energy that redshifts as Erad ~ This is 
the crucial feature that prohibits strings from dominating the energy density 
of the universe during the radiation era and makes the cosmic string theory 
a viable theory for structure formation in the universe. 

3.2 Scaling Solution 

Using the Monte-Carlo initial conditions described above, the Nambu action 
to evolve string segments and intercommuting to describe intersection events, 
it is straightforward (though not easy) to construct numerical simulations de- 
scribing the evolution of the string network in an expanding universe. Using 



this approach it has been shown |70| that the initial string network quickly 
relaxes in a robust way to a scale invariant configuration known as the scaling 
solution. According to the scaling solution the only scale that characterizes 
the string network is the horizon scale at any given time. On scales larger 
than the causal horizon scale t, the network consists of a random walk of 
long strings which are coherent on approximately horizon scales. On scales 
smaller than the horizon more recent simulations |[7l|, [72| , |73| (Fig. 7) have 
shown that the network consists of a fixed number of approximately 10 long 
strings coherent on horizon scales and a large number of tiny loops with typi- 
cal radius 10~^t. The efficient formation of tiny loops by the intercommuting 
of long strings on small scales leads to the existence of wiggles on the long 
strings (wiggly strings). The effective mass per unit length of wiggly strings 
is larger than the bare mass (/ie// — 1.4/i) and their tension is smaller than 
the bare tension. The main features of wiggly strings |7^, [75| , [76| , |77 
discussed in more detail below. 



are 
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Figure 7: The string network in the matter era. The volume of the cube 
is and the simulation box was plotted after an expansion by a factor 

a = 16 (from Ref. |7T[). 



There is a simple heuristic way to understand why does the network of 
strings approach a scaling solution with a fixed number of long strings per 
horizon scale. Consider for example the case when the number of long strings 
per horizon increases drastically. This will inevitably lead to more efficient 
intercommuting and loop formation thus transfering energy from long strings 
to loops and reducing the number of long strings per horizon back to its 
equilibrium value. Similarly if the number of long strings decreases reduced 
intercommuting and loop formation will tend to increase the number of long 
strings towards an equilibrium value. These heuristic arguments have been 
put in more detailed form using differential equations in Ref. ||29|| . 

3.3 Benefits of Scaling 

The scaling behavior described above is a crucial feature of cosmic strings. It 
ensures that the string energy density remains a fixed fraction of the matter 
density and thus an overclosure and premature collapse of the universe is 
avoided. To see this we may calculate the total energy in long strings at a 
time t assuming a fixed number of M long strings per horizon volume. The 
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energy density of strings is 



Pstr ^ ~ (^PPc (34) 

where I used the definition of the critical density pc ~ (^)^ (see Eq. (3) with 
A; = 0). Eq. (34) imphes that 

~i G/x ~ 10-6 (35) 



P 



c 



for rj = lO^^GeV. Therefore the density in strings remains a small fixed 
fraction of the energy density of the universe. As will be seen, this nice feature 
is not shared by domain walls even if it is assumed that walls approach a 
scaling solution. 

Gauge monopoles and domain walls are inconsistent with standard cos- 
mology each for a different reason. In the case of gauge monopoles there is 
no long range interaction and therefore there is no efficient mechanism to 
convert energy from monopoles into radiation. The monopole-antimonopole 
anihilation for example is very inefficient without long range interactions. 
Therefore, during the radiation era monopole energy redshifts with radiation 
as pmon ~ fl-^ and scales as a{t) leading to Pmon{to) » Pcito) and a 
premature collapse of the universe [[7^]. 

If Grand Unification to a simple group is realized in nature, the formation 
of monopoles is inevitable. This is the well known monopole problem of 
standard cosmology and may be seen as follows: Consider the symmetry 
breaking G — > of a simple GUT simple group G to a group H. In any 
realistic theory, H must include the gauge group of electromagnetism U{l)em 
since it is an experimental fact that the photon is massless and therefore 
f/(l)em is unbroken. Now, the vacuum manifold at the broken symmetry 
phase is M = G/H. The homotopy sequence of homotopy theory may be 
used to show that 

7r2{G/H) = 7T,{H) = 7,,{U{l)em) = Z (36) 

and therefore since tt2{M) ^ 1 monopoles must form in GUTs. 

The dilution of the monopole density during an epoch of exponential ex- 
pansion of the universe (infiation) provides one solution to the monopole 
problem]^. An alternatively solution is provided by constructing GUT 



21 



models where monopoles get temporarily connected by U{l)em strings during 
a temporary breaking of electromagnetism. Such models have the disadvan- 
tage of being somewhat unatural by introducing extra phase transitions in 
the symmetry breaking sequence of GUTs. 

GUT domain walls are inconsistent with standard cosmology |pO|| even if 



they manage to achieve a scaling solution. Consider for example a domain 
wall network with a single wall spanning each horizon scale. The wall energy 
density may be easily found as 

Pu> (37) 

Thus for t = to and rj ^ lO^^GeV we obtain 

^(to) ^ 10'' (38) 

Pc 

The wall symmetry breaking scale must be at least 17 orders of magnitude 
smaller than the GUT scale for a model containing walls to be viable without 
inflation. 

The above discussion shows the unique attractive features of cosmic strings 
compared to other defects in the context of standard cosmology. Several other 
such features will be discussed in the rest of this review. 



4 Gravitational Effects 

4.1 String Metric: The Deficit Angle 

The most important interaction on cosmological scales is gravity. It is there- 
fore important to understand the gravitational effects of strings before at- 
tempting to study in more detail their cosmological effects. The straight 
string solution is thin, cylindrically symmetric and Lorentz invariant for 
boosts along the length of the string. This imposes the following constraint 
on components of the energy momentum tensor T^^, 

T:{p) = T°(p) ^ f,6{x)6{y) (39) 

Also 

= ^T:{p) = O^T^ = Ty = (40) 
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where use was made of the cyhndrical symmetry and of the fact that 

^ as r — > oo. Therefore, the string energy momentum tensor may be 
approximated by 

T^^ ~ S{x)S{y)diag{ij, 0,0,-11) = diag{p,p^,Py,p^) (41) 

which imphes that the string has significant negative pressure (tension) along 
the z direction i.e. 

Pz = -p = -fJ' Px=Py = (42) 

This form of T^jy may be used to obtain the Newtonian hmit for gravitational 
interactions of strings with matter. For the Newtonian potential \Ef we have 

V^^ = 47tG{p + Y,Pi) = 0^ Fn = (43) 

i 

and a test particle would feel no force by a nearby motionless straight string. 

Simulations have shown however that realistic strings are neither straight 
nor motionless. Instead they have small scale wiggles and move with typical 
velocities of Vg — 0.15c coherent on horizon scales. What is the energy 
momentum tensor and metric of such wiggly strings? 

The main effect of wiggles on strings is to destroy Lorentz invariance along 
the string axis, to reduce the effective tension and to increase proportionaly 
the effective mass per unit length of the string. Thus, the energy momentum 
tensor of a wiggly string is 

T^, ~ 6{x)6iy)diag{peff, 0, 0, -T) (44) 

with T = —pz < p, Peff > p and PeffT = p^ {p is the 'bare' mass per 
unit length obtained from the field Lagrangian). As expected the above T^^/ 
reduces to the straight string case for p^fj = T. The breaking of Lorentz 
invariance by the wiggles also induces a non-zero Newtonian force between 
the wiggly string at rest and a test particle since the tension (negative pres- 
sure) is not able to completely cancel the effects of the energy density (Eq. 
(43)). This may be seen more clearly by using the Einstein's equations with 
the tensor of Eq. (44) to find the metric around a wiggly string. The result 



in the weak field limit (small Gp) isffT? 



ds' = (1 + hoo){dt' - dz' - (1 - ^GpeffYr'd^') (45) 
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with 

hoo = 4G(/ie// - T)ln{r/ro) (46) 

where tq is an integration constant. Clearly, in the presence of wiggles (/Xe// 7^ 
T) the Newtonian potential hoo is non-zero. The change of the azimouthal 
variable (f to the new variable ip' = (1 — 4G'/ie//)v9 makes the metric (45) very 
similar to the Minkowski metric with the crucial difference of the presence 
of the Newtonian hoo term and the fact that the new azimouthal variable (p' 
does not vary between and 27r but between and 271 — SirGfieff- Therefore 
there is a deficit angle^^^ a = SnG^eff in the space around a wiggly or 
non-wiggly string. Such a spacetime is called conical (Fig. 8) and leads to 
several interesting cosmological effects especially for moving long strings. 



Figure 8: The gravitational effects of a long string: The string deficit 
angle a leads to sharp discontinuities in the temperature of the CMB, velocity 
fluctuations in matter and formation of wakes, and lensing of galaxies and 
quasars. The Doppler and Sachs- Wolfe effects (see section 6.2) due to plasma 
velocities and potential fluctuations on the last scattering surface (tree) are 
also illustrated. 
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4.2 Velocity Fluctuations 

The main mechanism by which strings create perturbations that could lead 
to large scale structure formation is based on velocity perturbations created 
by moving long strings. Long, approximatelly straight strings moving with 
velocity Vs induce velocity perturbations to the surrounding matter directed 
towards the surface swept in space by the string. This effect may be seen 



more clearly by using cartesian coordinates!^ . 

Consider a straight long straight long string moving on the y — z plane 
with velocity Vs (Fig. 8). By transforming the line element of Eq. (45) to 
'shifted' cartesian coordinates defined as 

y' = r sin*/?' ~ y — r cosy? 4GyUe//V2 (47) 

x = r cos if' c::ir cos if = X (48) 

we obtain 

ds^ = {l + /loo) {dt^ - dz^ - dx^ - dy'"^) (49) 

where y' = y — AGfieff^px. The geodesic equations for a test particle in the 
cosmic string spacetime are of the form 

2x = -(I - ±^ - y'^)d,hoo (50) 

2y' = -{l-x'-y'^)dyhoo (51) 

Using the form of /iqo from Eq. (46) and perturbing around the initial particle 
trajectory (on the string frame) 

X = Vgt y' = Vo = const (52) 

the geodesic equations become 



X = -^(l-t;,2)4G(/ie//-T)^ (53) 
y = -^(l-t;D4G(/ie//-T)^ (54) 



A heuristic solution to these equations may be obtained as follows [ffTf : The 
Newtonian force induced by the wiggly string on a particle at distance r has 
a magnitude (cf Eq. (53-54)) 



25 



and is effective for approximatelly 



At ~ — (56) 
Tfierefore tfie magnitude of the induced velocity will be approximatelly 

and by symmetry it is directed towards the surface swept by the string. The 
exact solution to the system of Eq. (53-54) turns out to be very similar 

y'^y- 4.G;..,,..7. - '""^^^^^^ " (58) 

and therefore the total velocity perturbation induced by a moving long string 
to surrounding matter close to the surface swept by the string is 

A^; = ^^^^^-^^ ~ + iTTG^leffVsls (59) 

where the first term is due to the Newtonian interaction of the string wiggles 
with matter while the second term is an outcome of the conical nature of the 
spacetime. 



4.3 Microwave Background Fluctuations 

Another particularly interesting effect of the string indTiccd deficit angle is 
the creation of a characteristic signature on CMB photons. The type of this 
signature may be seen in a heuristic way as follows: Consider a straight long 
string moving with velocity Vg between the surface of last scattering occuring 
at tree and an observer at the present time to (Fig. 8). The metric around 
the string is approximated by Eq. (49) and at large distances from the string 
core we have 

dy' dy — AGneff^pdx — dy — AGfieff<^Vslsdt = dy — V{(fi)dt (60) 

For photons we have ds'^ = and therefore the Newtonian term hoo has no 
effect. Thus the presense of the moving long string between the observer 
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and the last scattering surface induces an effective Doppler shift to the CMB 
photons. For photons reaching the observer through the 'back' ('front') of 
the string we have 

(^)i,(2) = (^)i,(2) = ±{V{v = 27T) - Vi^ = vr)) = ±AnGfi,ffVs^, (61) 

where the 1(2) and +(— ) refer to photons passing through the 'back' ('front') 
of the string. 

Therefore a moving long string present between tree and today induces 
line step-like discontinuities on the CMB sky with magnitude [p2|, p3| 

Ti - T2 

— - — = SnGneffVs'ys (62) 

Notice that there is no Newtonian term inversely proportional to the wiggly 
string velocity as was the case for the induced velocity perturbations. In 
section 6 it will be shown that Eq. (61) can be used to construct a superpo- 
sition of long string perturbations thus constructing approximations to the 
predicted CMB maps. 

4.4 Lensing 

The existence of the deficit angle in the cosmic string space-time implies that 
strings act as gravitational lenses with certain characteristic properties 



p3| , 0. A long string present between a quasar and an observer (Fig. 8) 
will lead to the formation of double quasar images for the observer. The 
separation angle of the two images depends on the magnitude of the deficit 
angle and may be found by simple geometrical considerations. Using Fig. 8 
the separation angle AO is obtained as 

(/ + d) sin(^) = sin(|) ^ ^ SvrG/ie//^ ^ 5" (63) 

Therefore a gravitational lensing event induced by a cosmic string is expected 
to involve a number of neighbouring double images with typical separation 
of a few arcsec. Such a candidate event has indeed been observed and it will 
be discussed in some detail in section 5.5. 
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5 Structure Formation 



5.1 Sheets of Galaxies, Filaments 

In the cosmic string model for structure formation, small and intermediate 
scale structure (galaxies and clusters) is seeded by string loops while large 
scale structure is produced by perturbations induced by long strings. 

A loop with radius R has a mass Mi = f3Rn where /3 is a parameter 
approximatelly equal to 27r and on distances much larger than its radius it 
produces a gravitational field which is identical to the field of a point mass 



with mass M;[£5|. Therefore string loops can act as seedlike perturbations 
leading to the formation of clusters and galaxies. Recent simulations [[7T|, [72[ 
however have shown that the typical size of loops is probably too small to 
have a significant effect on the formation of objects like galaxies or clus- 
ters. This implies that the simple 'one loop one object' assumption that was 
made in the early days of the cosmic string model is probably incorrect and 
more sophisticated methods must be developed to study the way small and 
intermediate structure forms in the string model. 

The indicated relatively small importance of loops further amplifies the 
cosmological role of long strings. These strings can lead to large scale struc- 
ture formation through the velocity perturbations, produced during their 
motion, to surrounding matter ||8^. As discussed in section 4 long strings 
moving with velocity Vg produce velocity perturbations directed towards the 
surface they sweep in space. These perturbations whose initial magnitude is 
given by Eq. (59), grow gravitationally and within approximatelly a Hubble 
time t they form planar density perturbations called wakes |8^, |86| , |88| , pO| 
(Fig. 8). The gravitational growth of wakes can be calculated analytically 
in the linear regime using a simple but powerful method known as the the 
Zeldovich approximation. Using this method the thickness and typical di- 
mensions of the dominant predicted planar structures can be calculated |P9|| 
and the result can be compared with observations of redshift surveys in order 
to test the cosmic string model. In what follows I will sketch the basic steps 
of the calculation involving the Zeldovich approximation. 

Consider a long straight string moving with velocity Vg and sweeping 
a plane in an expanding universe. Consider also a test particle located a 
physical distance h{t) from the plane swept by the string. The scale h{t) 
will initially grow with the universe expansion but due to the string induced 



28 



velocity perturbation the growth will be decelerated by gravity, the scale h{t) 
will stop expanding, it will turn around and collapse on the string induced 
wake. In order for the planar structure formed to be able to fragment and 
lead to the formation of galaxies and clusters within it, it is necessary that 
the developed overdensity be nonlinear i.e. ^ ~ 0. It may be shown that 
this condition is indeed realized within the scales hni{t) that have turned 
around at a given time t. Thus, the scale hniit) that turns around at t is 
called the nonlinear scale at time t and the time when a scale h turns around 
to collapse is called the time of non-linearity tni{h) for the scale h. Clearly, 
hni{t) defines the thickness of the planar structure at the time t. The length 
of this structure is approximated by the coherence length of the string which 
is given by the horizon scale t while its width is approximated by the distance 
Vst travelled by the coherent portion of the string within an expansion time 
t. 

Observations indicate that the typical thickness of the observed sheets 
of galaxies is about 5h^^Mpc. By demanding hni{to,Gfi) = 5h^^Mpc i.e. 
that the predicted thickness of string induced sheets of galaxies is equal 
to the observed, the single free parameter of the model G/j, can be fixed. 
This normalization can then be compared with others coming from other 
observations and also from microphysical constraints. The length and width 
of the dominant structures can also be obtained. 

In the context of the Zeldovich approximation the physical (Eulerian) 
coordinate r of a test particle with initial comoving position q (Lagrangian 
coordinate) is written as 

f{q,t)=a{t){q-^{q,t)) (64) 

— * 

where ^ is the comoving displacement induced by the initial perturbation. 
In order to find the comoving scale qni{t) that turns around at time t we 
must first use dynamics to find "if and then solve r = to find qniit). This 
calculation is outlined below. 

The dynamics of the Eulerian coordinate r may be obtained in the New- 
tonian approximation using the equations 

'Kit) = -^W.t) (65) 
— $(f,t) = 47rGp(f,t) (66) 
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In addition, mass conservation implies that the mass in a Eulerian volume 
d^r should be equal to the mass in the corresponding Lagrangian volume 
a^(Pq i.e. 

p{r,t)dh^a{tfpo{t)d\ (67) 

or 

p{r,t) = a(i)Vo(t)^^ ^ PoW(l + ^^^m) (68) 
Prom Eq. (66) and Eq. (68) we obtain to hnear order in ^: 

^$(f, t) ~ 47rG[^po(i)r + po{t)a{t)^ {q, t)] (69) 

Using the Friedman equation Eq. (3) in Eq. (69) and using also Eq. (64-65) 
we obtain 

^ + 2-^ + 3-^ = (70) 

a a 

which with the appropriate initial conditions can determine the evolution of 

The initial conditions corresponding to the velocity perturbation induced 
by a moving long string at an initial time ti are of the form 

^{ti) = (71) 

^xiti) = ^y{U) = Q (72) 

A^; = a{ti)^^{U)^AT:GpVs-fsf (73) 

with / = 1 + 2{.l^^/y implied by Eq. (59). 

The growing mode solution of Eq. (70) with initial conditions (71-73) in 
the matter era (a ~ i^/^) can easily be found to be of the form 

m - \^vC-)njf'% (74) 

It is now straightforward to find the comoving scale qni{t) that becomes 
nonlinear and turns around at time t 

j^{a{qni - ^{t)) = ^ qni{t) = 2^{t) (75) 
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Using also the observational fact that qni{to) ~ 5h ^Mpc we may obtain a 
constraint on the only free parameter of the model 

Gfieff ~ 0.7{vs^sf)-' X 10-^ (76) 

where ti = t^q has been used in Eq. (74) since fluctuations can not grow 
before t^q due to pressure effects present in the radiation component. This 
result compares favourably with constraints coming from a completely dif- 
ferent direction: microphysics. In order for GUTs to be consistent with 
low energy experiments and with constraints on baryon lifetime, the scale of 
GUT symmetry breaking rj must be of the order lO^^GeV. For strings formed 
during a GUT phase transition this implies that 

{G^l)GUT ^ {Gv')gut ^ 10-' (77) 

which compares well with the corresponding constraint from macrophysics 
(Eq. (76)). This striking agreement, which was first observed in studies 
of galaxy formation by string loops|91|, between constraints coming from 



completely different directions is an exciting example of a meeting point 
between cosmology and particle physics and has also been one of the most 
attractive features of the cosmic string model. 

The earliest string wakes forming at tgq are the most numerous and also 
they have the longest time to grow. Therefore they give rise to the dominant 
sheetlike structures by today. Their typical dimensions are determined by 
the size of the comoving horizon at t^g and by the nonlinear scale qni{teq,to) 
which determines their present thickness (Eq. (75)). Thus the predicted 
dimensions of these structures are 

et:r x ^^.t^r x g„K^e„ to) ^ C40 x v^AO x SM^c^ (78) 



where ^ ~ 1 is the coherence length of long strings. The above predicted 
dimensions compare reasonably well with observations assuming relativistic 
strings. Simulations have shown that a large portion of long strings have 
relatively small velocities Vg — 0.15 on horizon scales. Such strings will tend 
to form structures that are more filamentary than sheetlike. 

The above discussion on the form and dimensions of the predicted large 
scale structure has not taken into account the effects of free streaming and 
therefore it is valid only for CDM. The effects of free streaming present in the 
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case of HDM can be taken into account by preventing the growth of scales 
smaller than the comoving free streaming scale 

Af = v,{t)tla{t) (79) 

where the velocity of HDM particles {eg neutrinos) starts droping like 
l/a(t) after T ~ my when they become nonrelativistic. A sketch of the time 
dependence of A™*" is shown in Fig. 9. 



Figure 9: The time dependence of the free streaming scale \^°^. 

An important value for Aj is its maximum value A™""^' which for adia- 
batic fluctuations produced during inflation determines the minimum scale 
that can form independent of larger scales (without fragmentation). Adi- 
abatic (but not seedlike) perturbations on all scales smaller than A™"^ are 
erased by free streaming and structures on these scales can only form by 
fragmentation of larger objects. This is a problem in adiabatic perturbations 
with HDM because observations indicate [^] that galactic and cluster scales 
formed earlier than larger scales. This is difficult to explain in models where 
fragmentation is the only mechanism to form smaller scale structures. Seeds 
like cosmic strings survive free streaming and therefore smaller scale fluctu- 
ations in models with seeds + HDM are not erased but their growth is only 
delayed by free streaming Thus galaxies and clusters can in principle 



form independently of large scale structure in these models. For a neutrino 
with mass rriy = 25eV (enough to produce Q = 1 for h = 1/2) we find 

XJ^'- ~ A,(te,) ^ Gh^o'Mpc (80) 
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and therefore the growth of all scales less than Gh^^Mpc is delayed by free 
streaming. Thus the introduction of HDM in the string model has two main 
effects. First it delays the growth of smaller scales thus transferring more 
power on large scales and second the delayed growth tends to increase the 
required value of Gfi for nonlinear structures to form by today. A detailed 
calculation shows ||89[| that 

{GfJ.)HDM ^ 2 X 10-^ > iGfl)cDM (81) 



5.2 Peculiar Velocities 

Using the result of Eq. (74) for the comoving displacement it is straight- 
forward to find the magnitude and coherence length of peculiar velocities 
u{t(),ti) produced by long strings at an initial time U as observed at the 



present time to. The typical magnitude of these velocities isp5 



u{to, U) = ^{to, U) = ^At;(^)i/3 {U > teg) (82) 

where Av = AirGfieffVs'jsf ■ The coherence length of these velocity fields is 
given by the coherence length of the long string that produced them which 
in turn is about equal to the comoving horizon scale L(ti) = tf^ = J^^ ^ ~ 

1/3 

. Thus the magnitude of the predicted velocity field with coherence scale 
L is ^ 

u{to, L) c::^ 300^Qh—p-km/sec (83) 

where /xg = G/i/10~^, L^q = tg°''" = 14/i~^Mpc. This result is not in good 
agreement with peculiar velocity observations on large scales which indicate 
the presence of velocity fields with magnitude of about 600km/ sec coherent 
on all scales from 10h~^Mpc to larger than 50h~^Mpc. Even though it is 
possible to play with the normalization of the model to induce agreement with 
observations on a particular scale, the scaling does not allow agreement 
on a large range of scales. This problem of the cosmic string model which also 
appears in the standard adiabatic CDM model can be addressed by assuming 
velocity bias which essentially means that the observed velocity fields do 
not represent accurately the underlying velocity fields of dark matter since 
luminous matter evolves differently than dark matter. The motivation for 
such a conjecture however is not particularly appealing especially on these 
large scales. 
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5.3 Power Spectrum 



An important quantity that characterizes a pattern of fluctuations —{x) is 



p 



the power spectrum. Consider a Fourier expansion of a fluctuation pattern 

— {x)= I d^k^{k)e'^^ (84) 
p J p 

— * 

The power spectrum P{k) of the pattern is then defined as the ensemble 
average 

P(fc) =< > (85) 

It may be shown that the correlation function C{x) of the pattern smoothed 
on a scale Iq ~ /cq ^ can be obtained as the Fourier transform of the power 
spectrum i.e. 

C{x) =< ^(fi)^(fi + x) >^i= / d^kP{k)e'^-^W{k - ko) (86) 
p p J 

where W{k — k^) is a filter function which filters out the scales that do not 
contribute due to smoothing and finite volume effects. Thus, the contribu- 
tion of scales smaller than Iq ~ /cq ^ to the variance of the fiuctuations is 
approximated by 

< {^{x)f C(0),„ = 5(^0, tf klP{ko) (87) 

Now the existence of the scaling solution implies that string induced per- 
turbations Ss have a fixed magnitude on horizon scales at any given time. 
Thus 

^si^k = ^, ^) = constant (88) 

where Afc = o,{t)^ is the physical smoothing scale. At a later time t in the 
matter era the fiuctuations will have grown to 

5ik,t)^i^f/'5iXk^ti,t,) (89) 

But 

t- = Afc - a{ti)/k ~ k'^ (90) 
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and therefore combining Eq. (88) and Eq. (89) we obtain 6{k, t) ~ k"^ which 
combined with Eq. (87) gives 



P{k) ~ k 



(91) 



known as the scale invariant i9t] Harisson-Zeldovich power spectrum. This 
is a generic type of spectrum which is also predicted by models based on 
inflation!^. Any large variation from the scale invariant spectrum is cosmo- 
logically unacceptable since it would lead to ultraviolet or infrared disasters 
in the magnitude of fluctuations. 



Figure 10: The log of the power spectrum vs the log of the scale for 
strings (dashed) and gaussian scale invariant fluctuations. In each case the 
curve with more power on small scales corresponds to CDM, the other to 
HDM. (From Ref. H). 

Given that growth of fluctuations starts at ieg, perturbations on scales 
larger than the horizon at t(.q enter the horizon after t^.^ and they keep growing 
without delay thus keeping their power law form P{k) ~ k. On the other 
hand fluctuations on smaller scales enter the horizon at earlier times and 
their growth is delayed until teq when matter dominates. The delay is larger 
the earlier the scale enters the horizon i. e. is larger for smaller scales. Thus, 
we expect a bending of the spectrum at smaller scales |^8| to a power low 
P{k) ~ /c" with n < 1. An additional bending is expected to occur in the 
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case of HDM|]99| due to the effects of free streaming which delay (erase) 
seed (adiabatic) fluctuations on small scales. In fact in the case of adiabatic 
fluctuations we have not only a bending but a cutoff of the spectrum at the 
maximum free streaming scale. The above discussion is illustrated in Fig. 10 
where the spectra of adiabatic and string perturbations are plotted for CDM 
and HDM models. The attractive feature of the model of strings with HDM 
is the transfer of power to large scales without completely erasing the power 
on small scales. This is in agreement with observations showing more power 
on large scales than predicted by the standard adiabatic+CDM model and 
also allows for the independent growth and formation of structures on smaller 
scales. Detailed N-body simulations of the string+HDM model attempting 
to explore in more detail these features are currently in progress. 



5.4 Gravitational Radiation 

A particularly interesting cosmological constraint that can be imposed on 
the cosmic string model is based on the gravitational radiation produced by 
oscillating massive loops |1T00| , |101|| . These oscillations are expected to lead 



to a stochastic gravitational wave background which could be detectable as 
perturbations of the period of msec pulsars [ p,02|| . The phase ip{t) of the 



periodic signal detected by a msec pulsar can be expanded as 

^{t) =^o + ^t- ]^(pe + ipR{t) (92) 

the first three terms can be modeled based on observations and pulsar physics 
while ^nit) is called the residual phase and is due to perturbations of the 
signal period induced by either gravitational waves or by other sources (pulsar 
intrinsic noise etc.). The quantity that characterizes the energy density of 
gravitational waves of angular frequency uo is 

. (93) 

Po 

where pg is the energy density in gravity waves. For the stochastic back- 
ground produced by string loops it may be shown that 

p Ott 
fi,(u;, Gp) = 2.5 X 10-V6(^)' < ^UT) > h-\-f (94) 
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where 



1 /■*+^ 

< y^UT) >=< - ^Idt > 



(95) 



and P is the period of the signal in sec. The presently observed time residual 
is P < (fj^ > /27r ~ 1.5/isec and this implies an upper bound for /xg whose 
precise value depends on the string simulations. The most recent careful 
study of constraints from the msec pulsar [p!0T| gives a bound on /xg of /ig ~ 
2 X lO^^h^^/^. This result is based on the numerical simulations of Allen 
and Shellard[^] and it is projected that if the residual timing noise remains 
constant the bound in 1998 will be /ig ~ 5 x 10~^h~^^^ thus effectively ruling 
out the model. 



5.5 A String Detection? 

The detection of cosmic string segments through the detection of multiple 
lensing pairs of galaxies or quasars is a very exciting prospect. The expected 
properties of such lensing events can be summarized as follows 

• A number of ~ 3 of galactic or quasar twin pairs is expected con- 
centrated in a small region in the sky [eg 50" x 50") 

• The angular separation of pairs is expected to be approximately con- 
stant of 0(1") as shown in Eq. (63). 

• In contrast to most typical lensing events strings are expected to induce 
no magnification of the images during lensing due to the conical nature 
of their metric. 

• The members of the lensed pair are expected to have very similar prop- 
erties (redshifts, spectra etc.) since they originate from the same object. 

• The redshift of the pair is expected to be relatively large to increase the 
probability for a string being present along the line of sight {z ~ 0.1). 

An event that effectively fulfills all the above properties was detected in 
1987 [|103|| by Cowie and Hu who detected 4 twin pairs of galaxies in a 50" x50" 



angular region in the sky with typical angular separation 3", no magnification 
and very similar properties. These properties are shown in Table 1 (from Ref. 
[|103|| ) while a contour plot of the twin pairs is shown in Fig. 11. 
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Table 1: Galaxy Properties (Ref. [M ^ 



Figure 11: A candidate string detection: Four galactic twin pair in a 
sky area of 50 x 50 arcsec (Ref. ||103|| ). 



The strength of the candidates however was reduced by a later publica- 
tion which showed that the images of the pair members do not match in a 
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satisfactory way when images are studied in radio band flO^ . Even though 
the issue is far from being resolved there is currently no clear evidence that 
these events are induced by a string lensing. 



6 Cosmic Strings and the Microwave Back- 
ground 

6.1 Angular Spectrum Basics 

A pattern of CMB temperature fluctuations is characterized by two classes of 
properties: properties of the angular power spectrum and statistical properties 
(probability distribution function etc.). I will give an introduction of the an- 
gular spectrum basics [ p^ ] and discuss the relevant predictions of the cosmic 



string model obtained by using a simple analytical approximation ||106| , |107 
The predicted statistical properties of the CMB fluctuations will be discussed 
also. 

Consider a photon scattered for last time on the last scattering surface 
at tree and reaching an observer with observational resolution at time 
to- One of the main sources of CMB fluctuations is the Sachs-Wolfe effect 
which is the result of gravitational potential fluctuations on the last scattering 
surface. Thus, an initially isotropic photon wavefront has to climb out of a 
potential 6"$ whose depth varies with direction. The resulting fluctuations 
(^)a, smoothed on a scale A on the last scattering surface are 

(^)a ~ S^x (96) 
The gravitational potential due to a mass overdensity 6p on a scale A is 

A 

where use was made of 6px ~ {k^P{k)y^'^ ~ ^(3+n)/2 ^ ^-{3+n)/2^ Now, Eq. 
(96) and Eq. (97) imply that 

(^)0 ~ S<^e ~ e(i-")/^ (98) 
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and therefore the correlation function at zero lag (A^ = 0) for angular reso- 
lution is 

C{A9 = 0)e i^-^r (^)L. - e(^-") (99) 

where <> denotes ensemble average (equivalent to angular averaging with 
the crgodic hypothesis). Therefore for a scale invariant power spectrum (n — 
1), {^)rms is independent of the angular resolution 

{—)rms ~ constant (100) 

Note that even for scale invariant primordial spectra (100) will be violated 
on angular scales less than the angular scale of the horizon at recombination 
{Qrec — ^reJ'^ — 0.03rad ~ 2°) due to microphysical processes taking place on 
subhorizon scales (Doppler effect etc.). Such processes are discussed below. 

CMB fluctuations are usually measured on large parts of the sky and 
in some cases over the whole celestial sphere. Therefore a convenient ba- 
sis to analyse these perturbations is not the Fourier basis but the spherical 
harmonics. A fluctuation pattern ^{q) can be expanded as 

^(g) = E«rw0) (101) 

l,m 

Using the addition theorem, the lack of preferred direction and defining 
Ci =< |a™p > we may expand the angular correlation function C{9) us- 
ing the angular spectrum coefficients Ci 

m ^< ^(gi)^(g2) 2^E(2/ + mPii^ose) (102) 

where cos^ = qi ■ q2. For observational experiments taking place on a rela- 
tively small part of the sky the contribution of low values of I is filtered out 
and the above expression reduces as expected to a two dimensional Fourier 
transform with the angular spectrum Q playing the role of the power spec- 
trum P{k) 

-,00 

c(^) = :r E {2i + i)CiPiicose) ^ — E iCiMW) ^ (103) 

C{9)eo ^ J dHCi/'W{l - lo) (104) 
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where W{1 — /q) is a filter function centered on tlie maximum sensitivity 
mode Iq and filtering out modes that are undetectable due to either limited 
sky coverage or limited resolution. 

Consider the power law ansatz Q ~ The angular scale corresponding 
to the mode / may be approximated by 9/ = j. Therefore, using Eq. (104) 
we can express {^)rms in terms of the resolution 60 and comparing with Eq. 
(99) we can express a in terms of the power spectrum index n i.e. 

C(0)eo = (f^)Ls - llQo - ©o"-' (105) 
Now comparing Eq. (99) with Eq. (105) we obtain a = n — 3 and therefore 

Q ~ (106) 
which is approximately valid for / > 1. Fig. 12 (from Ref. ||115| ]) shows 



the form of the angular spectrum Z^C; as predicted by some models based on 
adiabatic fluctuations. 



Figure 12: The CMB angular spectrum as predicted by some typical 
models based on gaussian scale invariant fluctuations. Using appropriate 
filter functions for each experiment the predicted (^)rms can be obtained. 

Notice the flatness of the spectrum {I'^Ci = const) on angular scales less 
than the horizon at recombination (Z ~ 100) which indicates the presence 
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of scale invariance. COBE, with angular resolution Oq — 10° and all sky 
coverage is sensitive to / ~ 35 and therefore the COBE data can provide 
both the normalization and the primordial spectral index n. Even though 
the different ways to perform the analysis of the data lead to slight variations 
in the value of n = 1, the data are consistent with a scale invariant spectrum 
n = 1 and seem to favor a small positive tilt (blue spectrum, n > 1) over a 
negative tilt (red spectrum tt, < 1). Most models based on inflation tend to 



favor a small negative tilt of the spectrum ||108| . 

Several experiments collect temperature data along one dimension in the 
sky for example along a meridian. The power spectrum index can also be de- 
rived by using the data from such one dimensional experiments. In order to 
avoid confusion with the two dimensional angular spectrum C;, I will denote 
the one dimensional angular power spectrum by P{k) and the variable con- 
jugate to the angle along the geodesic circle under consideration by k. This 
should not be confused with the density spectrum where a similar notation 
is usually used. The correlation function for a one dimensional pattern is the 
Fourier transform of P{k) 

C^\M) ~ — V P{k)e'^^'^W{k - ko) (107) 

where W{k — k^) is a filter function corresponding to the resolution of the 
experiment. By isotropy C^'^(A6') = C^'^(A^) and also k ~ I since they are 
both conjugate of the angular scale Therefore from Eq. (104) and Eq. 
(107) we obtain 

fCi ^ kP{k) ~ i^-^Us = C'io) = C'io) (108) 

These results will be applied in approximating the predicted angular power 
spectrum of the string model. 

6.2 String Angular Spectrum 

There are three main sources that can lead to CMB temperature fluctuations 
in the context of the cosmic string model: 



Kaiser-Stebbins [p2|, perturbations due to moving long strings present 



between tree and the present. 
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Potential fluctuations on the last scattering surface {Sachs- Wolfe effect |T09|) 



induced by long string wakes and loops with their gravitationally ac- 
creted matter. 

• Doppler perturbations induced by the local peculiar velocities of the 
plasma on which photons scatter for last time (see eg Ref. | [11(J| )). 

The above three types of perturbations are shown schematically in Fig. 8. 
Therefore the total string induced CMB temperature fluctuation may be 
wtitten as 

{—)tot = ijr)Ks + ijr)sw + ijr)D (109) 

Each one of these three contributions involves the superposition of some type 
of seeds and therefore in order to calculate it we must address the following 
two questions 

1. How can seeds be superposed? 

2. What type of seed should be superposed in each case? 

I will here briefly address the first question. A detailed study of both ques- 
tions may be found in Ref. [|1 1 1[] . For simplicity I will focus on one dimen- 



sional data. Consider a geodesic circle in the sky and a seed CMB tempera- 
ture fluctuation function fi{9) with amplitude ai and angular scale \l/ (Fig. 
13) centered at a random angular position 6i. 



Figure 13: Superposition of a seed on a geodesic circle on the sky. 
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A random superposition of N such seed functions will lead to a temper- 
ature pattern of the form 

N N 1 +00 

m = E «n/f - ^n) = E ^r.- E fl {k)e''^'-'-^ (110) 
n=l n=l A;=-oo 

and therefore the Fourier transform f{k) of the pattern may be written as 

N 

m^fnk)j:<^ne'''- (111) 

n=l 

Thus, the power spectrum corresponding to the pattern is 

Po{k) =< \f{k)\' >^N< |a„|2 > \f^{k)\' (112) 

In realistic cases however when the seed perturbations are induced by topo- 
logical defects, the size \& of the seed function will not be fixed but will grow 
as the horizon scale (the characteristic coherence scale of scaling defects) 
grows. Thus defect perturbations produced at later times when the horizon 
scale is larger will have a larger characteristic angular scale in the sky. As- 
sume for example that the comoving horizon grows by a scale factor a i.e. 
j.com _^ Q/^coTO ^j-^jg implies that the angular scale of the horizon Qh and the 
seed size ^ will grow by the same factor while the total number of seeds 
superposed along a circle at this later time will be smaller by the same factor 

^ aQh ^ (* ^ a*, N ^ N/a) (113) 

The power spectrum Pqik) corresponding to the resulting pattern after Q 
such expansion steps may be written as 

^Q(^) = E^i/r(^)r<Kr> im 

where the total number of steps Q may be obtained from the ratio of the 
maximum over minimum seed (or horizon) size while the total number of 
seeds at the first expansion step is the number of defects per horizon scale 
(obtained from the simulations leading to the scaling solution) times the total 
number of horizons present on the circle during the first expansion step 

27r 

-^max = a^^mm N = M- (115) 
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Eq. (114) may now be used to find the contribution to the total spectrum 
from each one of the following perturbation types 

• Kaiser-Stebbins (KS) perturbations 

• Potential perturbations due to wakes (W) and loops (L) with their 
accreted matter present on the last scattering surface. 

• Doppler (D) perturbations due to plasma velocities induced by strings 
(photons scattered on moving plasma suffer Doppler shift). 

The quantities that need to be specified in Eq. (114) for each one of the 
above perturbation types are, the type of seed function fi, the number of 
seeds per horizon M and the range of seed scales that need to be superposed. 
Assuming independent contributions from each perturbation type, the total 
spectrum Ptot may be obtained as a sum of the partial spectra as 

Ptotik) = Pxsik) + Pw{k) + Piik) + Poik) (116) 
The detailed derivation of the form of Ptot niay be found in Ref . ||1 1 1|| . There 



it is shown that Ptot{k) depends on four parameters; the single free parameter 
of the cosmic string model Gn which may be normalized by comparing with 
the CMB COBE data or by large scale structure observations and three other 
parameters which may in principle be fixed by comparing with numerical 
simulations. These parameters are defined as 

b = M<{vsisf> (117) 



/ = 1+ ^ (118) 

e : vl> = ^e,(t)/2 (119) 

where ^ is the string curvature radius \E' as a fraction of the horizon scale 
Qh{t) and / determines the wiggliness of long strings as defined previously. 
These three parameters may be fixed either by comparing directly with string 
simulations or by comparing the CMB spectral contribution P^sik) of Eq. 
(116) with the corresponding spectrum derived by propagating a photon 
wavefront through a simulated string network ||112| , |113|| p|. Both of these ap- 



proaches have been pursued in Ref. [ 111 | with results that are consistent 
with each other. 



^Ptot can not be used because CMB simulations with strings have not included so far 
the effects of potential and Doppler fluctuations 
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Figure 14: The total spectrum (a) including Kaiser-Stebbins, Sachs- 
Wolfe and Doppler fluctuations obtained as discussed in the text. The con- 
tribution of each individual component is also shown (b) (Ref. 



The total spectrum Ptot{k) in units of (G/i)^, obtained after the above 
normalization, is shown in Fig. 14a while the contribution of each type of 
perturbation is shown in Fig. 14b 

Doppler CMB fluctuations are due mainly to long strings present on the 
last scattering surface and therefore their characteristic scale corresponds to 
the coherence scale of those strings. By the scaling solution, this scale is 
of the order of the horizon at Uec {,Qh{trec) — 2° ^ khiUec) — 100). This 
explains the existence of a well defined peak for the Doppler term. The fact 
that the magnitude of this peak is significantly larger than the magnitude of 
the scale invariant Kaiser-Stebbins (KS) term may be understood as follows: 
The magnitude of the contribution to the KS term by each long string is 

ST 

{^)ks = 4:7TGnVs'jsk ■ {vs X s) (120) 

where k is the unit photon wave-vector and s is the unit vector along the 
string. The corresponding contribution to the Doppler term is 

(^)i, = k-6v = AnGfiVsjsk ■ {vs X s)f = {^)Ksf (121) 
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where / ~ 6 (/ defined in Eq. (118)) according to simulations ||71||. Thus we 
expect the Doppler term to dominate over the KS term for k ~ kh{trec) — 100 
as is in fact seen in Fig. 14b The neghgible contribution of loops (dotted line 
in Fig. 14b) to the total spectrum may also be understood by considering 
the fact that the typical loop radius is a tiny fraction (about 10~^) of the 
horizon as shown by simulations. Therefore a typical loop present on the last 
scattering surface at tree corresponds to an angular scale of about 0.3 arcsec 
which is way above the resolution of any present experiment. 

It is now easy to normalize the remaining parameter Gfieff by using the 
COBE DMR data. The angular correlation function predicted for the COBE 
(DMR) experiment can be expressed in terms of Ptot{k) and the window func- 

tion W{k — fco) — e (^-is)^ shown in Eq. (107). By demanding agreement 
with the detected {^)rms we have 

(^)fif = (C(0))^/^ = 1.1 X 10-^ {G^\ff = 1.6 X 10-6 (122) 



Figure 15: The cosmic string predicted correlation function smoothed 
on COBE scales. Superimposed are the first year COBE data (Ref. ||111|| ). 

This result is consistent with previous analytical studies IT 06|| and with 
numerical simulation studies ||112| , |113|| . Fig. 15 shows C{9) obtained by Eq. 
(107) (normalized with Eq. (122)) superposed with the COBE data. As 
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expected by the scale invariant nature of the string perturbations on COBE 
scales, the agreement is fairly good. 

The predicted power spectrum index can be found with a best fit based 
on Eq. (108). By using modes with k < 20 and the spectrum of Fig. 14a we 
find 

(123) 



n ~ 1.35 ±0.5 



Eq. (123) is a lower limit to the actually predicted n because Eq. (114) gives 
a slight overestimation of power on large scalesQ by assuming a fixed at all 
expansion steps (the expression 0h(t) — z{t)~^^'^ was used for all redshifts 
z{t) even though this approximation starts breaking down at low redshifts). 

Table 2:Detections of ^^^^^ x 10® and the corresponding predictions of 
the string [Qq = 1, h = 0.5, no reionization, A = 0) and inflationary models 
(0.8 <n< 1.0, A = 0) normalized on COBE. 



Experiment 


ko 


Ak 


Detection 


Strings 


Inflation 


COBE 





18 


11 ± 2 


11 ± 3 


11 ± 2 


TEN 


20 


16 


< 17 


13 ± 3 


9 ± 1 


SP91 


80 


70 


11 ± 5 


20 ± 5 


12 ± 2 


SK 


85 


60 


14 ± 5 


19± 4 


12 ± 3 


MAX 


180 


130 


< 30 ifiPeg) 


21 ± 5.5 


16 ± 5 


MAX 


180 


130 


49 ± 8 (GUM) 


21 ± 5.5 


16 ± 5 


MSAM 


300 


200 


16 ± 4 


19 ± 4 


24 ± 6 


OVR022 


600 


350 




13 ± 4 


17± 7 


WD 


550 


400 


< 12 


17.5 ± 4.5 


7 ±2 


OVRO 


2000 


1400 


< 24 


13.5 ± 3.5 


7 ± 3 



If the exact Qhi^) relation was used, a would need to be larger at late 
times thus reducing the number of expansion steps at large scales. The result 
would be a slightly reduced power on large scales and a tilt of the spectrum 
towards n > 1 for the scale invariant KS contribution. This effect was taken 
into account in Ref. [ p.l2|| where it was shown that it can increase the KS 
contribution to n by about 40%. Here, the effect on the total spectrum will 
be smaller since the other contributions remain unaffected by this. Using 
W{k) = e~^''~''°^ and fixing fco? for some of the ongoing experiments 



we are in position to predict the corresponding value of —^.^^ thus testing 
the cosmic string model. These predictions with la errors coming from the 



thank Neil Turok for pointing this out 



48 



variance of shown in Table 2. I also show some of the detections and 



114 


as 


U 


115 



models for 0.8 < n < 1.0, A = | |115|| ). At this time both inflationary models 
and cosmic strings appear to be consistent with detections at the la level. 
However, as the quality of observations improves, this may very well change 
in the near future. 



6.3 Statistical Tests 

The predictions of the string model on the CMB angular spectrum are useful 
in testing the model by comparing the amplitude of the predicted CMB 
fluctuations at various angular scales with the corresponding observations. 
However, as shown in Table 1, these predictions can not distinguish the string 
model from models based on inflation. A potentially interesting way that 
can lead to this distinction between models is the study of the non-gaussian 
character of string induced fluctuations. Before discussing a simple model 
that can lead to the identification of this non-gaussian character I will give 
a brief introduction to some statistics basics [|116|] that will be used later. 

A pattern of fluctuations of a random variable 6 is characterized by the 
probability distribution P{6) that gives the probability that a value 6 will 
be detected after a random sampling of the pattern. The n*'* moment of the 
distribution P{S) is defined as 

< 5" >= J d5 5''P{5) (124) 

The most useful moments are =< 5 > (the mean), cr^ =< 6"^ > (the vari- 
ance), as = (the skewness) and 04 = (the kurtosis). A distribution 
is completely defined by its infinite set of moments. 

The moment generating function is a very useful function that can pro- 
duce by differentiation all the moments of a distribution. It is defined as 

Ms{t) = Jd6 e'^P{6) (125) 

From this definition it immediatelly follows that 

<'^">=17;7^WI*=o (126) 
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An interesting example is the gaussian distribution defined as 



Pi6) = ^=e~^ P(5') = ^e-^ (127) 
V27rcr^ V2vr 



where 6' = ^— ^ is a standardized random variable. The generating function 
corresponding to the standardized gaussian distribution is obtained from Eq. 
(127) and Eq. (125) 

M5/(t) = e*'/2 ^ ag = 0, a4 = 3 (128) 

A particularly useful theorem that will be used below states that the moment 
generating function of a sum of two independent random variables is the 
product of the generating functions corresponding to each variable 

Ms,+s, = Ms.Ms, (129) 

Another basic theorem of statistics is the Central Limit Theorem (here- 
after CLT) which states that the sum 6n = xi + ... + x„ of a large number 
of independent random variables Xi with identical distributions P{xi) has a 
probability distribution that approaches the gaussian in the limit of n — oo. 

The above concepts apply in an interesting way to the CMB temperature 
fluctuations. The primordial fluctuations produced during inflation are due 
to the quantum fluctuations of an almost free scalar field. Therefore they may 
be represented as a Fourier series where the phase of each mode is random 
and uncorrelated with others i. e. 

^ k 

where 6k is a random phase. Therefore, since 5 is a sum of an infinite number 
of random variables, by the CLT it will obey gaussian statistics. 

The CMB fluctuation pattern induced by cosmic strings may be expressed 
as a superposition of seed functions as in Eq. (110). In what follows I will 
focus on the KS term and find the probability distribution P{S) of fluctua- 
tions induced by this term [|107| , |118|, |119| , |120|] . The result obtained in this 



way will only be an upper bound to the non-gaussian character of fluctua- 
tions because by the CLT any additional random effect can only make the 
fluctuations more gaussian. 
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Consider a step function perturbation with spatial size 2\1/ superposed at 
a random position on a pixel lattice with size 2L with periodic boundary 
conditions. The probability that the temperature of a pixel will be shifted 
positively (negatively) after this superposition is p = * (p = Thus the 
moment generating function of the fluctuation probability distribution after 
this single superposition is 

M,,{t) =< e*^i >= (pe* +pe-* + (1 - 2p)) (131) 

The corresponding result after the superposition of n seed functions at ran- 
dom positions is 

M5„^,,+,..+,„ = (2pcosht + (1 - 2p)Y (132) 

where use of the theorem (129) was made. Assuming small t (since to get 
the moments we set t = after differentiation) and p while taking large n we 
may write 

M,„ ^ (1 + ^^)" - e^'/' (133) 

where o"^ = 2pn is the variance of the gaussian distribution obtained (as 
expected by the CLT) for a large number n of superposed seeds. 

The moments of the distribution may now be easily obtained from Eq. 
(132) with the help of Eq. (126) as 

= ^Msm^o^^np (134) 

as = —Ms{t)\t=o = (135) 

a. = ^^^(^)l- = 3+1^^3 (136) 

and the gaussian value of the kurtosis is approached as 1/n. Clearly the 
crucial quantity that determines how close we are to the gaussian probability 
distribution is the product np. By setting P = j; and n = (number 
of seeds per horizon times the number of horizons in the lattice) we obtain 
np = M ~ 10. With this value, the relative deviation of the kurtosis from 
the gaussian is 

«4-<— < 1/2M^, 1 ^^3^^ 



3 60 
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and any non-gaussian feature is undetectable. Can np be decreased in order 
to amplify the non-gaussian features of the pattern? Clearly, in order to 
change n we would have to change the model under study since n is deter- 
mined by the scaling solution whose parameters are fixed for cosmic strings. 
In order to change p we would only have to change the type of seed function. 
This is indeed possible by considering the probability distribution of the tem- 
perature of a linear (or non-linear) combination of neighboring pixels. The 
most convenient combination is the temperature difference of neighboring 
pixels since this leads (for step-like temperature seeds) to a very localized 
seed function of 5-function type. The number of affected pixels (and there- 
fore p) is minimized by such a seed function. Let us therefore calculate the 
moment generating function corresponding to the variable d = (5*+^— 5* where 
the superscript now denotes pixel location. By defining a new probability 
q = J that a pixel will be affected by a given seed we have (in analogy with 
Eq. (132)) 



e*^i >= (2ge* + 2qe~^ + ge^* + ge"^* + (1 - 6g)) (138) 

and following the same steps as above we can obtain the moments for the 
new variable d as 03 = (due to the symmetry of the seed function) and 

04 = 3 + 1^1^ (139) 
Anq 

But now q — j; and n — which implies 

M 

nqr^ — <^ M = np (140) 
and the relative deviation of the kurtosis from the gaussian is 

04 ~ 3 ~ 12M 6M ^ ^ 

The above described test based on the CMB temperature gradient is effective 
if the the number of pixels ^' affected by half the step function is much larger 
than 2 (see Eq. (141)). Therefore we require 

* = > » 2 (142) 
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where is the pixel size in degrees. Thus the statistics of the tempera- 



ture gradient patterns provide a more sensitive test for non-gaussianity than 
statistics of temperature patterns ||117|, [I18|j in experiments with resolution of 



56pj^ ~ 1° or better. 



A more detailed analytical study of the moments of the string induced 
CMB fluctuations taking into account the growth of the seed functions as 



the horizon expands may be found in Refs. ||118| , p.07|| Those analytical re- 
sults can be tested by using Monte Carlo simulations performed as follows: 
A one dimensional array of 512 initially unperturbed pixels is considered and 
a number = M-^^ of step-like seeds i'^min is the initial seed size in units 
of pixel size which is determined by the resolution of the simulated exper- 
iment) is superposed on the array at random positions and with periodic 
boundary conditions. To simulate the growth of the horizon a larger seed 
size \1/ = a'^rnin is considered next and the superosition is repeated for Q 
such steps until the seed size becomes much larger than the size of the array 
of pixels (larger seeds shift the whole lattice by a constant and do not affect 
the statistics). The resulting pattern is then Fourier transformed and gaus- 
sianized by assigning random phases to the Fourier modes and reconstructing 
the pattern. The gaussian and the stringy patterns are then compared using 
different statistical tests, various resolutions and adding random gaussian 
noise of various signal to noise ratios. The Monte Carlo pattern constructed 
along the above lines may be written as 

Q N/a" 

^ wi.tr = E E«"/r*(^-^.) (143) 

^ q=l i=l 

while the Fourier transform is obtained as 



T^^'^^Jl I-l ^^^^^ 



and the gaussianized pattern with the same scale invariant power spectrum 
is 

S.rp +00 

V^l^— = E \-^me''^e^'-'^/^ (145) 

fc=— oo 

where 9k is a random phase. 
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Figure 16: An exaggerated stringy temperture CMB pattern composed 
of temperature Kaiser- Stebbins discontinuities(a) compared with its random 
phase reahzation (b) (Ref. ||119[| ). 



Figure 17: Plot of kurtosis of temperature fluctuations versus angu- 
lar resolution {dmin) with zero noise. Even at the lowest resolutions (~ 
0.5arcmin), this test cannot distinguish the gaussian from the string pat- 
tern (Ref. [[TT9|] ). 
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Fig. 16a demonstrates a stringy pattern with unrealistically amplified 
non-gaussian features (small number of seeds, fixed seed size) wliile Fig. 16b 
is tlie corresponding gaussianised (random phase pattern). These two types 
of patterns may now be compared using statistical tests based on the mo- 
ments of the distributions. After constructing 100 realizations in each case, 
the mean and standard deviation (cosmic variance of string perturbations [1T2I| ) 
of the kurtosis is found for various experimental resolutions starting from 0.5 
arcmin to one degree. Even without noise added this test proves inadequate 
to identify the non-gaussian features of the stringy pattern (Fig. 17). 

This could have been expected in view of the above analytical considera- 
tions. The corresponding test for the temperature difference of neighboring 
pixels however is much more efficient as shown in Fig. 18a where the kurto- 
sis of the two types of patterns are compared. For resolutions better than a 
few arcminutes the non-gaussian features of the stringy pattern are manifest 
through a significant increase of the kurtosis. 



Figure 18: Plot of kurtosis of temperature gradient of fiuctuations versus 
angular resolution. For small resolutions and zero noise (a), this test can 
distinguish between seed and gaussian patterns. With gaussian instrumental 
noise added (b) noise to signal ratio of n/s = 0.3 destroys all trace of the 
non-gaussian features of the seed pattern. The signal survives at n/s ~ 0.1 
(Ref. Ilig). 
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Unfortunately this nice feature of this test is very sensitive to the addi- 
tion of noise to the fluctuation pattern as shown in Fig. 18b which shows 
the effect of adding a modest amount of noise (the noise to signal ratio was 
nj s = 0.3) on the kurtosis of the temperature difference pattern. This sen- 
sitivity to noise of the temperature difference patterns is due to the low 
temperature derivative for most pixels in temperature patterns constructed 
by step-function superposition compared to the much larger derivative of 
noise patterns. Thus it is much easier for noise to dominate in temperature 
difference (derivative) patterns. 

7 Conclusion 

The purpose of this review was to give an introduction to the basic features 
of the cosmic string model for large scale structure formation and to show 
some of the recent developments in testing the model by comparing with new 
detailed observational data. The main points stressed are the following: 

1. Cosmic strings are linear concentrations of trapped energy predicted 
by GUTs to form during a phase transition in the early universe. 

2. For a natural value of a single free parameter ^ (mass per unit length) 
strings may provide the seeds for structure formation. 

3. The main achievements of the model include the following 

• Three independent types of requirements lead to the same value 
for the free parameter ^ of the model: Large Scale Structure, 
Microwave Background, Grand Unified Theories. 

• There is a well defined mechanism for the formation of the ob- 
served sheets and filaments of galaxies on large scales. 

• The model can resurrect massive neutrinos as dark matter candi- 
dates since strings-|-HDM appears to be a viable model for struc- 
ture formation with amplified power on large scales (seeds survive 
free streaming). 

• The statistics and power spectrum of the predicted CMB fluctua- 
tions are consistent with current CMB observational data. 
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• The existence of a scaling solution makes the model consistent 
with standard cosmology in contrast with other models based on 
topological defects. 

4. The main challenges of the model include the following 

• The complexity of performing detailed simulations of structure 
formation including both the string network and matter. This 
results to a large uncertainty in the predictions of the model. 

• The scaling of the predicted velocity fields with scale L is 
hard to reconcile with large scale velocity field observations. 

• The msec pulsar constraints appear to uncomfortably close to rul- 
ing out the model. 

The model makes some well defined unique predictions that can be 
used to rule out or verify it 

• The CMB fiuctuations are predicted to have certain non-gaussian 

features on scales of 1' or smaller. 

• The model predicts the concentration of multiple images of lensed 
objects in small areas in the sky with certain unique characteris- 
tics. Such an event may have already been detected. 
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